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1 . INTRODUCTION 

1.1  Statement  of  the  Problem 

The  frequently  occurring  meteorological  anomalies  in 
certain  geographical  areas  of  the  world  can  modify  significantly  the 
electromagnetic  propagation.  The  modification  in  propagation  affects 
the  performance  of  the  communication  systems  in  such  a fashion  as 
to  be  able  to  communicate  for  distances  significantly  larger  than 
what  is  expected  during  normal  meteorological  conditions.  For 
ground  based  tactical  missile  systems  and  the  long  range  surveillance 
radars , the  impact  is  more  significant.  The  radar  systems  experience 
range  and  angular  errors  for  look  angles  which  are  greater  than 
a couple  of  degrees.  For  elevation  angles  less  than  a couple  of 
degrees,  the  systems  can  detect  targets  at  significantly  large 
distances.  However  the  signals  corresponding  to  these  long  ranges 
of  detection  are  contaminated  with  the  clutter,  the  primary  and 
multiple  time  around.  Further  along  with  the  extended  detection 

■ 

ranges,  for  slightly  higher  elevation  angles  the  target  detection 
is  impaired  due  to  the  existance  of  holes  or  the  shadow  zones. 

The  existence  of  regions  of  extended  detection  ranges, 
and  of  regions  of  shadow  zones,  are  dependent  on  the  meteorological 
conditions  which  tend  to  modify  the  refractive  index  profile  and 
hence  the  propagation  behavior.  The  refractive  index  rj  is  related 
to  the  atmospheric  pressure  P,  absolute  temperature  T,  and  the 

I 

partial  water  vapor  pressure  e.  The  refractive  index  n is  given  in 

terms  of  different  parameters  such  as  N,  N,  and  B units,  etc.  In 

I 

this  report  we  use  only  the  N and  M units  to  describe  the  refractive 
index  and  they  are  defined  as  follows: 


I 
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*1  k2e 

N » (H  -1)  x 106  * P + 

T <j2 

Here  the  constants  kl  and  k2  have  the  following  empirical  values. 

kl  - 77.6,  k2  - 3.73  x 105 

For  normal  atmosphere  N decreases  uniformly  with  height  imparting  a 
downward  bending  for  a ray  path.  On  the  other  hand,  the  earth's 
curvature  gives  the  appearance  that  the  rays  tilt  upwards.  This 
apparent  upward  tilt  can  be  nullified  by  incorporating  the  curvature 
of  the  earth  in  the  refractive  index  profile  which  leads  to  a 
straight  trajectory  for  a flat  earth.  The  media  dependent  bending 
of  the  ray  then  is  described  in  terms  of  modified  refractive  index 

a 

profile  n or  in  terms  of  M - units,  and  they  are  defined  as  follows: 

~ z 

n * n + — 
a 

and 

A z 

M » (n  -1)  x 106  a (n  + — - 1)  X 106 

a 

Here  z is  the  height  and  a is  the  radius  of  curvature  of  earth, 

(a  * 6.4  x 10*>  m)  . 


f 


Thus  M and  N units  are  related  to  each  other  in  the  following 
fashion: 


M - (n  + ( — x 106) ) 
a 

For  normal  meteorological  conditions,  the  plots  of  M and  N units 
as  a function  of  height  reveal  that  there  is  a steady  increase 
of  N units  with  height  as  contrast  to  steady  decrease  of  N units 
with  height  and  this  is  displayed  in  Fig.  (1-1). 
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During  anomalous  meteorological  conditions  the  temperature  and  the 
water  vapor  content  change  in  such  a fashion  as  to  cause  significant 
deviations  in  the  refractive  index  profiles  and  thus  in  the  profiles 
of  M and  N as  a function  of  height. 

Some  representative  plots  of  M profiles  have  been  re- 
produced in  Pig.  (1-2)  and  they,  in  general,  fall  in  three  major 
categories: 

1.  Normal  refraction 

2.  Subrefraction 

3.  Super refraction 


200  300  400  M-Units 

Figure  (1-2)  - M - Profiles  for  (1)  Normal,  (2)  Subref ractive 
and  (3)  Superref ractive  Conditions 

The  subrefraction  and  superrefraction  cases  are  due  to  the  anomalous 
meteorological  conditions  and  this  report  deals  with  the  electro- 
magnetic propagation  through  such  an  anomalous  medium. 
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1.2  The  Outline  of  the  Report  j 

This  •action,  Section  1,  has  identified  the  main  topic 
of  this  report,  namely,  the  determination  of  propagation  through 
meteorological  ducts. 

Por  the  propagation  studies,  the  meteorological  data  for 
Canton  Island  was  provided  by  ESD  and  Section  2 evaluates  the  data 
and  classifies  it  in  an  appropriate  fashion. 

The  electromagnetic  propagation  through  the  medium 
around  Canton  Island  can  be  studied  by  geometrical  optics  which 
provides  a reasonable  characterisation  of  the  line  of  sight  pro- 
pagation. Section  3 discusses  this  approach  and  its  limitations. 

The  alternate  approach,  namely  the  physical  optics.  Is  also  discussed 
in  this  section  in  order  to  highlight  its  application  to  the 
dlfraction  region.  Section  3 further  provides  the  ray  trajectories 
for  four  days  of  meteorological  conditions  so  chosen  as  to  cover 
the  broad  spectrum  of  variations  in  the  characteristic  parameters 
of  the  duct  such  as  the  duct  width  and  the  lapae  rate  of  the  N 
units. 

The  primary  purpose  of  this  investigation  is  to  develop 
the  physical  optics  approach  for  propagation  through  meteorological 
ducts  and  hence  details  of  the  analytical  formulation  of  the  guided  ‘ 
mode  approach  is  provided  in  Section  4 for  two  refractive  index 
profiles,  bilinear  and  trilinear.  The  numerical  analysis  based  on 
several  techniques  such  af  Newton- Raphaon  to  determine  the  characteristic 
wave  numbers  has  been  carried  out  for  the  bilinear  case  and  is 
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discussed  in  Section  4.  Also,  the  height-gain  functions  which 

describe  the  dependence  of  electromagnetic  field  with  height 

for  different  modes  have  been  included  in  Section  4 for  the  bilinear 


refractive  index  profile. 


Finally,  the  conclusions  and  recommendations,  the 


references,  and  the  different  appendixes  are  provided  in  Sections 
5,  6 and  7,  respectively. 
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2.  METEOROLOGICAL  DATA  FROM  CANTON  ISLAND 

For  the  purpose  of  anaysis,  it  was  decided  during  the 
early  stages  of  the  contract  that  Raytheon  will  be  provided  with 
sample  meteorological  data  from  Canton  Island  (28QS,  171. 6°W) . 

2.1  Data  Format 

The  data  was  provided  by  U.S.  Air  Forces  Environmental 
Technical  Applications  Centre  at  Scott  Aire  Force  Base,  Illinois  and 
was  in  terms  of  the  meteorological  parameters  (T,  P,  e).  Here  T,  P, 
e,  respectively,  stand  for  the  temperature,  the  atmospheric  pressure 
and  the  partial  water  vapor  pressure.  The  data  provides  for  the 
typical  refractive  index  profiles  for  the  year  1973.  This  data  is 
translated  into  N units,  M units  and  the  rate  of  change  of  N with 
height  in  tables  (2-1)  through  (2-12).  These  tables  are  for  twelve 
months,  with  two  profiles  per  month.  The  criterion  that  is  used  to 
determine  the  presence  or  absence  of  a duct  is  based  on  the  condi- 
tion that  the  rate  of  change  of  N units  with  height  is  less  than 
106/r  where  r is  the  radius  of  the  earth  (r  ■ 6.4  x 106m) . 

2. 2 Data  Evaluation  and  Classification 

Based  on  this  criterion,  the  following  observations  are 
made  regarding  the  nature  of  the  refractive  index  profile: 

1.  The  meteorological  duct  was  present  for  all  the  days 
for  which  the  data  has  been  provided. 
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32.200 

2861.64 

•5.00000 

16,325 

24,900 

2905.31 

•7.84983 

, , , t.  , - . - 

18.723 

27.200 

2965.06 

-6,78392 

20.730 

19,000 

3271.79 

•4.08570 

...  21.775 

15.900 

3432.66 

•2.96651 

1 1 — r - , — 

23.860 

10.700 

3754.62 

-2,49400 

24.530 

9,700 

3656.75 

•1.49254 

25.023 

8,600 

_ 3935.21 

•1.82556 

„ , n • — - - -| 

26.450 

7.100 

4157.42 

•1.19131 

26.776 

6.700 

4206.49 

•1.21951 

..  27.311 

6.100 

.4241,52 

•1.12570 

30.379 

3.600 

4770.63 

•,74967a 

j 

2MAR 

M 

. 0 

N 

- 29 

M 

ON/OH 

.003 

375,100 

375.571 

.062 

328,500 

338.229 

•769.831 

1.079 

258,100 

090.173 

•09,6820 

2.179 

229,600 

571.511 

•00.7103 

3.124 

201,600 

. 691.793 

•29.6296 

5.000 

107,800 

1067,30 

•19.6637 

. 6.601 

135.200 

1170,96 

•17.0000 

7.600 

120.900 

1313.03 

•10.3103 

6.357 

110,700 

1022,01 

•13.0702 

9,740 

45,300 

1623.62 

•11.1352 

11.020 

82.800 

1611.97 

•9.76563 

12.520 

69,600 

2030,10 

•8.80000 

..  10.330 

56,300 

« 2300.85 

•7.30807 

16.209 

03,600 

2593 , 26 

•6,61803 

16.600 

00,600 

2657,60 

•6.96056 

17.540 

30,500 

. 2740.58 

•6.70330 

16.679 

26.000 

2959.36 

•5.60107 

19.622 

22.300 

3132,61 

•5.33683 

1 20.600 

16,600 

3251.19 

-0.09871 

22.130 

10.100 

3066,56 

•3.07J90 

23.175 

12,000 

3608,03 

-2,00957 

23.770 

10.600 

3700,00 

•2,35290 

AR 

M 

0 

N 

15  - . . . 

H 

ON/hH 

. .003 

391.700 

..  392.171 

.075 

373.200 

380.966 

•256,900 

1.090 

26  6 , 0 0 

520.026 

-61.0517 

3.102 

220.000 

713.017 

-39.8058 

6.613 

1*5,000 

. 909.235 

•23.7933 

5.670 

1*5.800 

1076.87 

•23.5481 

7.610 

122.900 

1317.00 

-18.9080 

9.750 

99 . 900 

1620.79 

-13.0801 

11.050 

*1.800 

1815.66 

•10.0769 

11.625 

76,800 

1900.90 

•8.69565 

12.560 

69,000 

2039.81 

•8.30225 

10.390 

55.600 

2313,76 

•7.21311 

15.912 

06.100 

2502.88 

•6.37319 

16.757 

00,100 

2669, U7 

•7,10059 

16.929 

39.300 

2645.66 

•0,65116 

TABLE  (2-3)  -March  1973 


w ; 


17 


% 

SAP  R 

H 

0 

N 

27 

M 

ON/OH 

. .003 

390,900 

395.371 

.075 

309,000 

360.768 

•637,500 

0 

1.0*5 

267,000 

501.583 

•57.7465 

. . 3.105 

206,100 

699.586 

-36.9001 

0.307 

179,500 

. 655.320 

•22.8916 

• 

0.708 

167.100 

918.395 

•25.7796 

5.880 

108,500 

1071.14 

•17.0330 

7.610 

121.000 

1315.50 

•15.6647 

• 

0.150 

102,700 

1539.07 

•12.1114 

. 0,700 

05.500 

1623.82 

•12.2867 

11.010 

*3,000 

1811.00 

-9.52756 

0 

12.500 

69.000 

2031.30 

•9.06040 

10.310 

56,100 

2301.51 

-7.62431 

16.680 

40,400 

2657.69 

•6.62447 

0 

17.663 

33.800 

2805.34 

-6.71414 

....  17.871 

31.700 

2835.87 

-10,0962 

18.069 

29.100 

2927.11 

•4.34723 

0 

16.700 

27.500 

2961.75 

•6.92641 

10*810 

21.900 

3130.33 

-5,04505 

20.700 

18,800 

3266.86 

-3.48315 

0 

21.061 

17.100 

3324.96 

•4,46194 

23.950 

' * • 

10.400 

..  3768.04 

•2.335M 

I 


5APR  0 20 

M H H TABLE  (2-4)  - April,  1973 


.003 

3*3.600 

364.071 

* ...  «...  . . 



.085 

331,300 

344.638 

•637.805 

1.505 

251,200 

487.353 

•56.4085 

3.153 

201.900 

696,643 

•29.9150 

5.262 

149,000 

984,671 

•20.3414 

5.890 

147,600 

1071.81 

-18.1529 

7.630 

121.300 

• 1318.54 

•15.1149 

- 

9.760 

95, 7C0 

1627.16 

-12.01*8 

' 

9 

11.030 

83.200 

1813.94 

•9.84252 

12.530 

69,600 

2035.71 

•9,06667 

14.340 

59,100 

2306.22 

-7.45856 

1 

16,730 

3*. 800 

2664.94 

•6,82008 

| 

18.760 

28,000 

2971.67 

•5.81291 

18.925 

27,200 

2996,76 

•4,84846 

0 

20.3*0 

20.100 

3217.97 

•4,87973 

I 

20.730 

18.400 

3271.19 

-4,85714 

. j 

23.960 

10,500 

3770.11 

•2.44582 

0 

26.610 

6.900 

0182,33 

•1.35849 

27.898 

5,500 

43*3.03 

•1,08696  

— „ 

0 

29.515 

4,400 

4635.66 

-.680272 

0 

I A 


'.'I  ■■ 

1»  - ' » 


fMAY 

H 

0 

N 

16 

32 

M 

ON/DM 

• 003 

379.000 

379.971 

,0*0 

393,900 

356.561 

-377.778 

.6*7 

312.200 

951.361 

*36.6593 

1.1?6 

279,700 

951.363 

-156.909 

1.500 

256.200 

. 969,196 

•9J.6508 

2.1*6 

220.000 

563.010 

•56.0117 

3. 1S« 

197.200 

692.726 

•23.9568 

0.398 

197,600 

670.053 

•7.89916 

• .709 

170.700 

922.151 

•38.7755 

5.690 

192.500 

1076.71 

•16.5309 

7,620 

120.900 

1316.07 

•16.5599 

9.790 

95.300 

1623.62 

•11.6396 

11.020 

•3.300 

1612.97 

•9,37500 

12.500 

70.300 

2031.70 

•8,76376 

19.300 

56.600 

2300,99 

•7.61111 

15.563 

96.000 

2990,0? 

-6,60916 

16*660 

39.900 

2659,05 

•7,36377 

17.766 

33,000 

2623.69 

•6,12769 

18.700 

27.000 

2961.25 

..  -6.56955 

19.515 

23.300 

3065.90 

•9,53988 

19.819 

21.600 

3131.99 

•5,5921 i 

. 20.730 

16,900 

...  . 3271.19 

•3.51262 

|0*4Y 

H 

0 

N 

23 

N 

ON/DM 

.003 

372,900 

372.671 

.065 

396.100 

361.936 

-296.391 

.797 

333,300 

950,513 

•22.3565 

.972 

300.600 

053.11* 

-105.333 

1.500 

279,500 

519.66* 

-39.9621 

1.839 

252.900 

591,061 

•78.9661 

1.150 

209.900 

696.673 

-36.9907 

1.527 

1*9,100 

792.529 

-90.5836 

9.566 

168.100 

667.699 

•19,6300 

5.9oo 

195,600 

1071.5* 

-16.9711 

7.620 

120.900 

1316,07 

•19.7670 

9.790 

•95.700 

1629.0? 

-11.6509 

11.010 

63.300 

• 1610.90 

-9.76378 

12.990 

70,900 

2030.23 

-8,71622 

15.2*1 

99.500 

2997.27 

•7,96836 

16,659 

00,200 

2653.91 

-6.77399 

19.6*1 

22.000 

3110.16 

•6,01255 

20.600 

18.700 

3282.97 

•2.99906 

29.026 

10.900 

37»0.37 

•2.57285 

?6.6«7 

6.900 

9105,99 

•1.31037 

30,167 

3.900 

9737.96 

-.869553 

31.933 

3,300 

9935.51 

-.973939 

32.909 

2,600 

5165.63 

-.975867 

TABLE  (2-5)  - May 


1973 


19 


3JUN 

H 

0 

N 

15 

H 

ON/OH 

.*•00  3 

360.600 

365.071 

.083 

336,700 

309.720 

•308.750 

l.ll* 

296.100 

O70.900 

•39.3792 

1 • o°7 

298.800 

083.697 

•123.099 

1.810 

•237.800 

, 522.038 

•30,7003 

1.950 

266.200 

572,178 

208.82a 

3.136 

216.100 

708.076 

•01.9099 

3.231 

1«5.500 

702.083 

•220,000 

5.559 

151.300 

1023.57 

•18.9863 

. 5.850 

155.200 

1073,10 

13.0021 

6.282 

107,600 

1133.32 

•17.5926 

6.600 

100.000 

1175,62 

•23.899a 

7.580 

120.600 

• 1309,99 

-19.7959 

9.680 

96.000 

1610,91 

•11.7103 

9,799 

98,000 

1631.98 

•13.0050 

5JUN 

H 

* 

0 

N 

23 

H 

0N/DH  TABLE  (2-6)  - June, 

1973 

. . ,.003 

371.600 

372.071 

.071 

304,200 

355.301 

•002.901 

.789 

320.000 

003.800 

•33.7007 

1.086 

266,300 

099.071 

•77.0005 

1.99J 

220.200 

. 536.926 

•83,0375 

2.970 

221.200 

687.856 

•3.05810 

3.131 

207.700 

698,°91 

. -85.9873 

3.79a 

1*3.500 

778.820 

•36.5008 

a.602 

175,100 

903.085 

•9.90566 

. 5.212 

160,800 

978.625 

•25.0877 

_ 

5.096 

160.700 

1027,09 

13.732a 

5.850 

105.900 

1063, eu 

•53.1073 

6.105 

109,700 

• 1107.65 

10.9020 

6.368 

103.000 • 

1102.22 

•25.0753 

7.580 

122.000 

1311.79 

•16.9967 

8,207 

115.000 

1009.05 

-10.0908  . 



9.690 

96.900 

1617.38 

•12.8205 

10.960 

83.500 

1803.26 

•10.5512 

*1  1 

12.000 

70.300 

2022.2* 

•6.91892  . _ 

10.230 

57.000 

2289,86 

•7.03017 

10.767 

53.000 

2370.12 

•7,00879 

. 15,006 

47.600 

2071.26 

•7.952*7 

. i 

16,590 

00,000 

2603.57 

•6.29371 

j 

| M | | 1 1 l - . ^ . 

3JUL 

H 

0 

N 

H 

ON/OH 

..  .003 

353,300 

353.771 

.075 

328.600 

390.366 

•393.056 

1 .«^1 

25a , aoo 

*66.767 

-52.7738 

1.6*2 

?39 , <>oo 

503.526 

•73.6318 

l.Ofll 

225.300 

. 529.666 

•55.2129 

3.1?a 

1 °7 , 300 

667. 693 

•23,6686 

. . 3.95a 

1*5.100 

805.530 

•19,6988 

a. 637 

166.600 

69a. 201 

•27,0669 

5.660 

1 a 5 , 9 o o 

1065. ao 

•16.9256 

6.905 

131.000 

12ia,o8 

-19.2S8U 

7.570 

121.500 

1309.32 

•19.2857 

6.033 

115.100 

1375.57 

•13. *229 

. 6.605 

1*6,900 

1956.63 

•11.7133 

6.663 

10a,  700 

1998.55 

•13.309a 

9.670 

96,000 

1613.39 

•11.0596 

..  10.802 

ea,9oo 

1779.66 

-9.00565 

1 0 , 9a  o 

83,500 

1600.12 

•10,1999 

12.«20 

70,700 

2019,55 

•8,69865 

• 

• 

9JUL 

H 

0 

N 

31  . ... 

M 

OV/DH 

.003 

370.300 

370.771 

.067 

305,900 

356.913 

•381.250 

.666 

307,800 

997.138 

-96.9068 

1.979 

237.300 

969.J73 

•119.289 

2.070 

222,100 

. 596.908 

•25.7191 

2.563 

210.200 

612.365 

-29.1379 

3.122 

197. aoo 

687.279 

-22. *980 

5.190 

157.900 

963.926 

•19.6216 

5.660 

106.300 

1065,60 

•15.9167 

7.570 

121.200 

1309,02 

-10.6769 

6.799 

1*5.600 

1966.27 

•12.6932 

9.660 

96, aoo 

1615.31 

•10.9U27 

10.990 

89.100 

• 1600.72 

•9,76190 

12.910 

71.000 

2018.2* 

•8.91156 

19.200 

56,900 

226S.05 

-7.87709 

16.006 

' 99,300 

255S.53 

•6,97e79 

16.590 

39.900 

2635.22 

-8.23970 

17.077 

36,900 

2715,90 

•6.51769 

17.338 

33.300 

2753.89 

. -1  1.8770 

16.690 

26.200 

2951,09 

•5.95315 

20.001 

20.800 

3159.20 

-3.96767 

20.710 

18.000 

3267,65 

-3.99922 

23.197 

12.300 

J699.39 

•2.33890 

23.950 

10.500 

3768.59 

•2.20159 

21 


2AUG 

H 

0 

N 

32 

M 

DN/DH 

. .003 

371.600 

372.071 

.003 

307.100 

361.693 

•272.222 

.910 

303.100 

005.890 

•53.8556 

1.500 

267.200 

502.566 

•60.6U75 

1.570 

235 , <100 

081.752 

-050.286 

1.826 

228.700 

515.221 

-26. 1719 

2.635 

209.100 

622.563 

•20,2270 

3.133 

197.500 

689.105 

•23.2932 

3.337 

192.900  ' 

716.515 

•22.5090 

5.538 

152,500 

1021.06 

-18.3553 

5.850 

107. 0OC 

1060,90 

-17.6282 

7.560 

121 .200 

1307.05 

-15.0877 

9.660 

96.900 

1612.67 

•11.5710 

10.521 

88,200 

1739.07 

•10.1005 

10.920 

63.700 

1797.18 

•11.2782 

..1 2. 3°0 

70,900 

2015.00 

•8.70708 

13.595 

61.000 

2100,62 

-7.88362 

10.170 

56,800 

2280.20 

-8.00000 

15.531 

• 

07,100  .. 

. 2060,10 

-7.12711 

SAUG 

H. 

0 

N 

32 

H 

DN/DH 

. ....  .003 

360,700 

361.171 

.0*0 

335.300 

306.081 

-313.580 

.509 

321,000 

007.1 05 

-30.7527 

1.223 

278.700 

070.603 

-62.7596 

1.983 

256.200 

066.900 

•86.5385 

3.111 

196.100 

686.253 

-35.6*80 

3.620 

1 *6 , 800 

750.621 

-22,2000 

0.627 

167,800 

693.632 

•16,8679 

5.820 

106,800 

1060.03 

-17,6027 

7.530 

121.000 

1302.95 

-10.8538 

8.970 

103.300 

1511.03 

-12.5306 

9,630 

06.000 

1607,06 

-10.5183 

10.600 

86,800 

1750.07 

-9.89601 

10,8®0 

■3,700 

1792.07 

-10,6807 

12.360 

71.000 

2010.03 

-8,63906 

13.017 

62.500 

2167.79 

-8.00163 

H.150 

56.700 

2277.00 

-7.91269 

is.3ie 

08 , 600 

2052.1* 

-6.93093 

_ 16*235 

01. coo 

25*8.07 

-8,28790 

16.510 

39,100 

2629,72 

•6.9*909 

16.992 

36,000 

2702.65 

•5.60166 

17.810 

30.200 

2825.03 

•7.50258 

18.630 

26,200 

2909,07 

-0.90196 

. 22 


a 


• 

9SEP 

H 

0 

N 

42 

M 

DN/OM 

1 

.003 

391,500 

391.971 

,00a 

370,200 

389,950 

•239,066 

• 

.3«3 

355, H00 

909.621 

•57.2313 

' 

.707 

321.500 

*32,437 

-99.2308 

1.501 

2«3.000 

52*. 925 

•39.3909 

• 

1 

• 

1.905 

262.700 

561.617 

-75.9901 

. . 2.515 

223.700 

612.33* 

-63.9399 

A 

2.876 

210.800 

662.079 

-35.73*1 

, 

o 

3.  l?5 

219,800 

710.150 

36.19*6 

3.07J 

205,600 

750.555 

-90.80*6 

( 

3.9*6 

196. UOO 

221 .*51 

-17.9337 

9 

0.116 

199,600 

**5.650 

26.1538 

- 0,380 

176,500 

• 863.77* 

-88.2576 

1 

5.380 

199,600 

1003.79 

-16.9000 

f 

# 

5,650 

162,500 

10*9,05 

10.7*07 

« 

5.830 

153,000 

1067.80 

. -52.7778 

i 

6.057 

1<16, UOO 

1159.58 

•10.5263 

1 

• 

7.530 

122.300 

1303.85 

•22. *609 

1 

6,330 

111.600 

. _ 

-1918.68 

-•13.3750 

1 

6.090 

112,800 

1995.61 

7.31707 

1 

6.599 

1*9,500 

1958,79 

-3l.02e6 

1 

8.726 

109,300 

1976.51 

-1.57a*0 

6.920 

109.900 

1509,05 

•25.2577 

j 

• 

9.620 

96,600 

1606.09 

-11.1*29 

- 

1 - - - - 

. ..  . 0.561 

*7,300 

— 

. 1799,95 

. _ -9.88310 

— — — - — . . - _ „ 

j 

o 

* 

I 

I 

o 

. 

’ 

5SEP 

0 

17 

A 

H 

N 

M 

DN/OM 

TABLE  (2-9)  - Sept.,  197  • 

1 

.003 

3*6,700 

369.171 

— .. ......  . . — * _ 

.102 

353.900 

369.905 

•3S1.515 

• 

# 

.619 

301.900 

*30.911 

-72.529* 

1.151 

2*6.700 

067.306 

•99.7831 

... 

1.99? 

296.790 

560.053 

-97.2813 

• 

2.232 

292.700 

592.92* 

-17.0213 

2.o*o 

223.900 

613.66® 

•79.6032 

■ — - - - - _ . . 

3.139 

197.600 

690.197 

-90.1527 

• 

0,799 

169,800  ‘ 

917.821 

-19,7590 

5,809 

196,900 

1069.6* 

-17.0476 

_ > . . . 

6,060 

116.600 

1150,8* 

-16.7990 

• 

7.565 

121.300 

1308.39 

•13.8965 

• 

9.690 

06.500 

* 1616.9* 

-11.6796 



10,607 

*7.200 

17®1.57 

. -10.1418 

# 

10.950 

89,0*0 

1802.  19 

-9,329*5 

• 

. . _.12.“20 

70.800 

2019.65 

-8.97959 



10.165 

57.300 

2279.96 

-7.73639 

• 

• 

I 


1 

70CT 

H 

„ o 

H 

37 

M 

ON/OH 

.005 

365.000 

385.971 

.078 

357,000 

370.139 

-361.333 

1 

.609 

320.500 

430.181 

•60.225a 

l.ooa 

282,200 

453.862 

-96,0620 

i.oos 

267,200 

. 500.214 

-3.8,3632 

1 • 

2.269 

217,300 

572.540 

-64.0S65 

..._  5.118 

107,600 

. 686.852 

-23.0670 

5.715 

1«6,500 

770.926 

-1.84255 

t • 

a. 021 

176,800 

609.743 

-57.6431 

..  ..  5.890 

195,000 

1062.27 

-16,0869 

6. 200 

130,600 

1126,56 

-14.0000 

6.591 

139.000 

1161,26 

-18.7251 

7.560 

121.500 

• 1307.75 

-13.1501 

8.105 

113.000 

1385.36 

-13.9063 

* 

0.670 

06,000 

1613.34 

-11.4231 

.10.790 

65,400 

1772.05 

-0.82301 

10.090 

83.700 

1600,32 

-8.90052 

t 

12.910 

70.800 

2018.06 

-8.77551 

14.200 

57.000 

2285.15 

-7.70950  

15.152 

50.400 

2424.79 

-7.08155 

• 

15,458 

47,400 

2472.95 

-9.20245 

— — 

16.004 

44,200 

2555.42 

•5.86061  ..  ..  . 
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TABLE  42-10)  Oct.,  1973 

A 

H 

K 

M 

DN/DH 

.005 

375.100 

375.571 

« 08a 

353.700 

366.681 

-264.108 

• 

.640 

317.500 

410.336 

-64,0706 

1.488 

263.  PO 

516.585 

-41.0012 

2.134 

259,800 

. 594.650 

-36.06*1 

• 

2.262 

238,600 

593.735 

•164,062 

2.980 

212.300 

601 .4a2 

-121,560 

3.108 

210,500 

696.182 

-2,8662a 

• 

5.830 

146.300 

1061.10 

-23.5856 

7.550 

121.000 

130S.69 

-14.7093 

0.660 

.06,000 

1611.77 

-11.8483 

• 

10.6*6 

85.000 

1762.66 

-9,64405 

10.020 

63.600 

• 1707,06 

-9,82906 

* 

12.400 

70.500 

2016.21 

-8,85135 

• 

14.189 

57.100 

2282.74 

-7.51121 

14.806 

52,700 

2375.90 

-7.07395  _ . 

* 

16.530 

40,500 

2634.25 

•7.07657 

• 

16.762 

38.000 

2668.16 

-10.7759 

18.274 

28,600 

2896.01 

_ -6, 21693  
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5N0V 

0 

20 

H 

N 

M 

ON/OH 

. ,003 

371.900 

372.371 

,067 

397.100 

357.613 

-387.500 

,696 

317,000 

926.211 

-97.8537 

1.968 

279,700 

505.097 

-59.7927 

1.958 

239.700 

. 593.795 

-79.9681 

3,107 

1*7,900 

660.02*; 

•36,1698 

S,  830 

196.300 

1061.10 

•16.7661 

7.550 

121.900 

1306.5* 

-19.1660 

8.683 

109, <100 

19*8.25 

-13.1283 

9.650 

*6.900 

1610.60 

•10.9302  . . . 

10.727 

85,500 

1768.6* 

•10.1207 

10.910 

83.700 

1795.61 

-9,83607 

12.3*0 

70.600 

’ 2019.79 

-8.85135 

19.180 

56,900 

2281.91 

-7.65363 

16.530 

90.300 

2639.05 

-7.06383 

17.329 

39.300 

2752.69 

-7.55668  . . . . 

16.592 

27.000 

2999.31 

•5.75710 

18.902 

29.800 

2990.75 

•7.09677 

....  20>5®0 

18,800 

3299.62 

. -3.55950  

21.516 

15,700 

0 

3391.82 

-3.39773 

* 

, TABLE  ( 2 — 1 1 ) - Nov.,  1973 


6N0V 

H 

_ 0 

N 

17  ..  .. 

M 

ON/OH 

.003 

363.600 

369.071 

.067 

392,500 

353.013 

-329,688 

1.966 

269,600 

999 , 9<|  7 

-52.0393 

1.696 

296.600 

513.037 

-100.000 

2.929 

218,700 

590.839 

-38.1669 

3.095 

209,900 

690.593 

-20.7207 

...  . 5.810 

150.200 

...  1061.86 

-20.1973 

6.163 

193,960 

1110.95 

•19.2635 

7.530 

121.000 

1302.55 

-16.3862 

9.690 

96,600 

1609.23 

-11.5690 

10.715 

85.600 

1766.91 

•10,2326 

10.9«0 

89,000 

1799.39 

•8,69865 

12.380 

70.500 

• 2013.07 

-9.12162 

13.298 

63,900 

2150.02 

-7.73420 

19.170 

86.900 

2279,89 

•8,02752 

19.957 

59.900 

2322.88 

•6.96869 

16.S90 

90.100 

2635,92 

•6.66510 

i 


7 


I 


SOEC 

0 

27 

ft 

H 

N 

M 

0*J/DM 

• 003 

370,000 

374.671 

,064 

3*55,200 

369.381 

-232.009 

» ' 

1.4*9 

287,500 

521.162 

•68.1851 

2.0«2 

200,200 

572.660 

•71.8076 

3.1U 

206,800 

605.626 

•36.5969 

3.035 

101.700 

800.160 

•18.3022 

4.41* 

175. 7C0 

86«.°37 

-33.1263 

5.379 

150.100 

008.126 

-22.6766 

5.810 

109,900 

1061.56 

•9,76478 

8,865 

13". 600 

1211.80 

•16,5026 

% 

7.131 

126,700 

1265.66 

-2^. 6992 

7.520 

121.600 

1301,38 

-13.6267 

0.810 

97.600 

•.  1605.32 

-11,6833 

• 10.009 

87.000 

1735.32 

•10,6862 

ft 

10.870 

86,100 

1780.73 

-10,2626 

..  12.3«0 

71.300 

2007,50 

•8.70768 

13.283 

63.800 

2166,06 

•7,95336 

• 

14.110 

57.300 

2271,33 

•7.85973 

0 


TABLE  ( 2-12)  - Dec.,  1973 


«DEC 

M 

0 
- N 

33 

M 

DN/OH 

.003 

.075 

366,400 

351.000 

366.871 

363.668 

•201.389 

.056 

200.200 

449.20* 

-59.8184 

. 1.676 

272.  100 

503.702 

-52.1154 

1.872 

260,800  ' 

. 554.534 

•2*. 5354 

2.328 

231.000 

596.291 

-65.3509 

3.101 

216.800 

703.384 

•16.3700 

•.205 

182,800 

856.737 

•26.9757 

9,436 

1*3,000 

682.200 

1.24224 

9.080 

161.800 

443.222 

-40.4580 

5.7*0 

152.500  • 

1059.45 

•11.6250 

6.223 

145.500 

1121,96 

•15,6014 

6.338 

138,200' 

• 1132,71 

-63.6763 

7.900 

123. 700‘ 

1298.97 

•12.5468 

8.500 

108,600 

1442,35 

•14.9505 

0.580 

46.700 

1599.92 

-11.0165 

10.642 

88,300 

1726,77 

•9, 7«478 

10.840 

84,600 

1785.53 

•9.29648 

12.300 

71,200 

2001.22 

•9, 1 7808 

12.063 

66,100 

2100.15 

•7.69231 

14,080 

*7,200 

2266,5? 

-7.96777 

15.810 

45.200' 

2527.39 

•6,90052 

16.480 

90,100 

2613.46 

-8.77797 
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2.  Eighty  percent  of  the  time  there  existed  surface 
ducts. 

3.  Twenty  percent  of  the  time  the  duct  profiles  were 
more  complex  with  possible  existence  of  elevated 
ducts. 

4.  The  limited  number  of  data  points  (M-units  versus 
altitude)  force  to  characterize  the  surface  ducts  as 
bilinear . 

5.  The  duct  width,  under  the  circumstances,  is  inac- 
curately determined.  For  the  available  data,  the 
duct  width  is  most  likely  the  maximum  value  rather 
than  the  likely  value. 

6.  The  duct  intensity,  as  determined  by  the  rate  of 
change  of  N with  height,  varies  from  -201/km  to 
-789/km.  The  minimum  intensity  occurring  in  December 
and  the  maximum  in  March.  The  refractive  index 
profiles  based  on  the  radiosonde  data  are  not 
corrected  for  different  measurement  errors.  The 
errors  could  either  be  sensor-related  or  could 

be  the  timelag  in  temperature  and  humidity  measure- 
ments. Finally,  the  most  conspicuous  error  is  the 
lack  of  correlation  between  the  first  radiosonde 


data  point  at  some  altitude  which  is  different  than 
mean  sea  level  and  the  sea  level  measurement. 
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7.  The  duct  widths  for  the  surface  ducts,  defined 
as  the  height  at  which  the  minimum  for  M occurs, 
varies  from  40  meters  for  January  2nd  to  102  meters 
for  September  5th. 

8.  For  the  largest  duct  width  of  102  meters,  the 
minimum  frequency  that  is  trapped  is  approximately 
350  MHz,  and  the  minimum  for  the  smallest  duct 

of  40  meters,  it  is  1.43  GHz.  This  is  based  on 
a crude  relationship  between  the  maximum  wave 
length  trapped  and  the  duct  width. 


'max  s*  0.014d3/2 
feet 
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cm 
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3.  GEOMETRICAL  VERSUS  PHYSICAL  OPTICS  APPROACH 

The  electromagnetic  propagation  through  a medium  is 
governed  by  the  refractive  index  profile  of  the  medium;  and  it  can 
be  characterized  either  by  the  geometrical  optics  that  is  the  ray 
treatment,  or  the  physical  optics  that  is  the  wave  treatment.  The 
distinction  in  the  two  approaches  which  was  highlighted  for  optics 
originally,  is  equally  valid  for  electromagnetic  propagation. 

3.1  Geometrical  Optics  or  Ray  Trajectory  Approach 

The  geometrical  optics  approach  is  primarily  based  on  the 
Fermat's  principle  which  states  that  the  line  integral  of  the  re- 
fractive index  between  two  points  is  stationary.  This  leads  to  the 
well  know  Snell's  law  which  states  that  niri  cos  0i  ■ n2^2  cos  $2* 

Here,  n^  ig  the  refractive  index  at  ri  and  0i  is  the  angle  of 
inclination  of  the  ray  with  the  horizontal,  and  n2  the  refractive 
index  at  r2  and  ?2  is  the  angle  of  inclination  at  r2. 

Based  on  this,  a set  of  ray  trajectories  for  actual 
refractive  index  profiles  has  been  plotted  to  highlight  the  salient 
features  of  the  geometrical  optics  approach.  The  ray  trajectories 
are  for  a specific  location  for  four  different  days  in  the  month  of 
August  at  local  noon.  The  antenna  height  is  50  meters,  and  the 
elevation  angle  is  between  0 and  1 degree.  All  the  four  ray 
trajectories  reveal  most  of  the  features  of  the  ground  duct,  namely, 
the  trapping  of  the  rays  at  angles  less  than  1°,  the  shadow  regions  beyond 
40  km  range,  and  the  modest  refraction  at  or  above  1°  elevation. 


A 
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These  ray  trajectories  reveal  the  cases  when  the  j 

geometrical  optics  approach  is  not  valid.  It  is  well  established 
that  the  geometrical  optics  or  ray  treatment  is  valid  when  two 
conditions  are  satisfied  and  those  are: 

1 An  1 j A Jl/2 1 

— <<  t <<  1 

kn  n kn  jl/2 

Here  J is  the  Jacobian  and  for  detail  explanation  see  Ref.  (2), 

(Paes  52  - 58)  . 

The  first  condition  states  that  the  relative  variation  of 
the  refractive  index  per  wavelength  is  smaller  than  unity.  The 
second  condition  states  that  the  relative  density  of  the  rays  in  a 
bundle,  or  the  relative  spacing  between  the  rays,  is  less  than  a 
constant.  This  second  condition  is  quite  often  not  adequately 
emphasized  when  the  geometrical  optics  approach  is  pursued.  This 
second  condition  is  often  violated  when  the  duct  is  present  as  is 
evident  in  the  ray  trajectories  shown  in  Figures  (3-1)  thru  (3-4). 

Further,  the  ray  trajectores  indicate  that  there  exists  a complete 
blank  region  with  no  electromagnetic  energy  at  all,  which  in 
practice  is  not  the  case.  Lastly,  the  geometrical  optics  approach 
is  independent  of  frequency  and  polarization.  This  aspect  prevents 
the  geometrical  optics  approach  from  treating  any  scattering  pro- 
blems and  thus  the  clutter  which  is  dependent  on  frequency,  pola- 
rization and  roughness  of  ground  or  sea  surface.  For  both  com- 
munication and  radar  systems  these  considerations  are  of  great 
importance,  and  the  preferred  approach  is  that  based  on  physical 


Figure  ( 3-3)  Typical  Ray  Trajectories  for  August  30 
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optics.  The  physical  optics  approach  provides  important  information 
on  the  extent  of  the  shadow  region  and  the  depth  of  the  "null"  in 
the  shadow  zone.  Further,  it  provides  a means  to  estimate  the 
clutter  whether  it  be  primary,  multiple  time  around,  or  sidelobe. 

3.2  Physical  Optics  or  Guided  Mode  Approach 

The  physical  optics  or  wave  approach,  which  does  not  have 
the  limitations  identified  for  the  geometrical  optics  approach,  is 
most  suitable  for  the  case  when  the  meteorological  ducts  are  intense 
enough  to  trap  the  waves.  The  approach  can  provide  system  related 
parameters  such  as  extended  detection  ranges,  field  intensity  levels 
in  the  shadow  region,  and  clutter  primary,  and  multiple  time  around. 
Thus,  in  the  diffraction  zone  (Fig.  3-5) , the  physical  optics 
approach  is  preferred  to  the  geometrical  optics  approach  when  the 
estimates  of  the  radar  performance  degradation  due  to  the  meteo- 
rological ducts  are  required. 

The  physical  optics  approach  essentially  consists  of 
determining  the  solutions  of  the  second  order  vector  wave  equation 
subject  to  the  appropriate  boundary  conditions.  For  electromagnetic 
problems  the  fields  E and  H,  which  are  solutions  of  the  Maxwell 
equations  subject  to  the  boundary  conditions,  can  be  described  in 
terms  of  scalar  modes  if  the  vector  field  can  be  resolved  into  its 
scalar  components.  This  is  highlighted  for  the  rectangular 
coordinate  system  in  the  following  equation: 


OtFFR.AtfcTtOM  2oM6 


Figure  (3-5)  - Tropospheric  Duct  Propagation 

jkxi*x 

£ (r)  - E a^i  (r)  • £ ai  ^ (T)  e 


Here,  £ ( r ) is  the  total  vector  electromagnetic  field  and  is 
described  in  terms  of  the  summation  of  the  characteristic  modes, 

(T) ; and  T stands  for  the  transverse  coordinates  y and  z which 
are  orthogonal  to  x,  the  direction  of  propagation.  The  longitudinal 
wave  numbers  kx^  are  the  characteristic  values.  The  modes  (T) 
are  functions  of  the  transverse  coordinate,  and  they  form  a complete 
orthonormal  set.  The  amplitude  coefficient  a^  thus  are  determined 
from  the  source  conditions  with  the  help  of  the  orthogonality 
condition.  Thus  implicit  in  this  approach  are  the  characteristics 
of  the  source  generating  field  such  as  frequency,  polarization, 
antenna  height,  etc.  Application  of  the  boundary  condition  implies 
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the  statement  of  the  field  characteristics  (a)  at  the  ground,  (b) 
at  the  different  regions  of  the  M profile,  and  (c)  the  Sommerfeld 
radiation  condition  for  altitudes  significantly  higher  than  the  duct 
heights.  This  approach  is  eminently  suitable  for  transmission  paths 
which  are  beyond  the  line  of  sight  or  in  the  diffraction  zone  and 
also  equally  applicable  at  the  turning  points,  and  in  the  shadow 
zone.  Thus,  since  the  objective  of  the  study  is  to  generate  radar 
coverage  diagrams,  when  the  meteorological  ducts  are  present,  the 
only  viable  approach  is  the  physical  optics. 

3.3  Ray  Trajectories  For  Four  Refractive  Index  Profiles 
The  geometrical  optics  or  the  ray  trajectory  approach 
involves  the  use  of  Snell's  Law  to  identify  the  propagation  paths. 
Typical  ray  trajectories  for  four  days  of  refractive  index  profiles 
from  Canton  Island  are  shown  in  Pigs.  (3-6  to  3-9).  For  January  and 
April  meteorological  conditions,  the  existence  of  shadow  zones  are 
clearly  noticeable.  This  should  be  evident  from  the  strength  of  the 
duct  as  measured  by  dN/dz. 

For  the  September  and  December  meteorological  conditions, 
the  duct  has  been  weak  and  a weak  shadow  zone  is  evident.  It 
is  to  be  noticed  that  the  shadow  zones  of  January  and  April  will 
not  be  as  pronounced  as  is  displayed  if  finer  intervals  in  the 
trajectory  angles  are  chosen.  Based  on  these  trajectories,  one 
would  conclude  that  in  the  shadow  zone  there  will  be  a total 
absence  of  propagation  paths  and  hence  a total  absence  of  electro- 
magnetic power.  The  question  which  is  very  pertinent  for  system 
performance  evaluation  is  how  valid  is  this  interpretation; 


Fig.  (3-7)  Ray  Trajectories  For  Canton  Island  - April,  1973 
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and  is  there  a better  way  to  estimate  the  field  strength  in  the 
shadow  zone. 

In  conclusion,  it  has  been  established  that; 

a.  the  geometrical  optics  approach  is  not  valid  in 
the  diffraction  zone; 

b.  it  is  inapplicable  when  the  relative  spacial  vari- 
ation of  the  refractive  index  per  wavelength  is 
larger  than  one; 

c.  when  the  fractional  change  in  the  spacing 
between  the  neighboring  rays  in  a wavelength  is 
not  less  than  unity,  then  it  is  not  applicable. 
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4.  GUIDED  MODE  APPROACH  - PHYSICAL  OPTICS 

It  is  thus  evident  that  the  guided  mode  approach  is  the 
most  viable  approach  to  establish  radar  performance  in  the  diffrac- 
tion zone  and  in  the  presence  of  meteorological  ducts.  The  measure 
of  the  radar  system  performance  is  either  (a)  in  terms  of  radar 
coverage  diagrams  depicting  signal  strengths  or  path  loss  as  a 
function  of  ground  range  and  altitude,  or  (b)  the  signal  to  noise 
ratio  as  a function  of  slant  range.  Such  a performance  evaluation 
is  mandatory  for  the  judicious  hardware/software  modifications  to 
alleviate  or  minimize  the  deleterious  effects  of  propagation  degra- 
dation. Recognizing  this,  the  goal  of  this  phase  of  the  program  is 
to  determine  propagation  behavior  in  the  presence  of  bilinear  and 
trilinear  types  of  refractive  index  profiles.  The  important  mile- 
stone of  the  next  phase  of  the  program  is  to  establish  the  radar 
coverage  diagrams  for  different  meteorological  conditions  and  for 
different  systems  parameters. 

4.1  Propagation  Through  Stratified  Medium 

For  the  propagation  analysis,  it  was  assumed  that  the 
medium  is  stratified  in  the  vertical  z direction.  The  stratifi- 
cation of  the  medium  dictated  by  temperature,  pressure  and  partial 
water  vapor  pressure  variations  with  altitude  is  described  in  terms 
of  the  refractive  index  profile,  which  in  turn  is  represented  by  the 
N-profile  or  the  M-profile.  The  profiles  which  have  been  considered 
are  the  bilinear  and  trilinear.  This  is  based  on  the  evaluation  of 
the  Canton  Island  meteorological  data  for  1973. 
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Height  ( Meters ' 


The  M-profile  Cor  January,  April,  September  and  December  are 
produced  in  Pig.  (4-1)  and  they  will  be  considered  in  the  numerical 
analysis. 


1 - Jan.  3,  1973 


Fig.  (4-1)  - Canton  Island  Refractive  Index  Profiles  for  Four  Days 
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4.1.1  Propagation  Through  The  Meteorological  Duct 
Witn  Bilinear  M Profile 

The  electromagnetic  propagation  through  any  medium 
is  governed  by  Maxwell's  equations. 

VxH-ju>eoE»0 


VxE  + ja,yoH»0 

and  the  individual  components  of  the  electric  and  the  magnetic 
fields  are  the  solutions  of  a second  order  vector  wave  equation 

(V2  + k2)  * (r)  = 0 

Here,  Laplacian  72  stands  for  the  following: 


2 fl  3 3 1 32  32  1 

U 3r  3r  r2  302  3Z2 J 

f 1 3 - 3 1 3 3 

— (R2  — ) + “ ; — (sin  e — ] 

3R  3 R R2  sin  9 39  39 

♦ i.  — — — ) 

R2  sin  9 3d>2  J 


V.3X2  3y2  3Z2  J 

in  cylindrical,  spherical  and  rectangular  coordinate  systems 
respectively. 


The  wave  number  k - kQ  n,  kQ  is  the  free  space  wave  number,  kQ  = 
u)/c,  u)  is  the  angular  frequency  and  c is  the  velocity  of  light  in 
free  space,  n is  the  refractive  index  of  the  medium. 

The  vector  field  9 stands  for  the  E and  H field 

/ E (r)  \ 

$ -►  ( ) 

\ H (r ) / 

For  a point  source  above  ground,  such  as  the  dipole,  the  wave 
equation  is  more  suitably  described  in  terms  of  the  spherical 
coordinate  system.  For  the  homogeneous  atmosphere  the  mode  solution 
of  the  wave  equation  was  first  obtained  by  Watson  (18)  and  is 
given  by 


in*.  “jkxm*x 

Ey  * e3“t  z — — -e  Umi  (z)  um2  (*> 

m xl/2 

Here,  umi  (z)  and  um2  (z>  are  the  height  gain  functions  at 
the  transmitter  and  the  receiver,  respectively;  the  subscript 
1 stands  for  transmitter  and  2 for  the  receiver.  The  transmitter 
and  receiver  locations  have  been  left  unspecified.  For  a ground 
radar  system,  the  transmitter  position  is  fixed  at  h and  then  Ujnl  (*) 
= um  (h)  is  the  specified  antenna  characteristics. 
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are  the  amplitude  coefficients  of  the  modes;  x is  the 
horizontal  range.  Finally,  the  height  gain  function  um  (z)  is  the 
solution  of  the  one  dimensional  wave  equation. 


f- 

l dz2 


+ kQ2  - kx2  \ Ey  (z) 


-jkx. x 


The  x dependence  has  been  assumed  as  e 

The  z dependence  will  be  deleted  from  hereon. 


and  Ey  (z)  = um2(z) 


The  refractive  index  n is  assumed  to  be  stratified 
only  in  the  vertical  z direction.  For  the  chosen  four  days  of 
M profile 


z < d 


z > d 


d is  the  duct  width. 


The  scalar  wave  equation  for  Ey  is  thus  given  by 


^■^2  + k02  (1  + 2 M (z))  - kx2J  Ey 


Thus  for  the  M-profile,  the  wave  equation  takes  the  following 
form: 


M (2)  * Mq  - pz 


z < d 


Mmin  + <J  <z-d) 


{^2  + Bl-C1Zl 

1 7z2  * B2  + C2  * If  % * ° 


z > d 


z < d 


z > d 


Here  Bi  . |<02  (1  + 2 Mq)  - kx2 


Cl  » 2ko2  P 


b2  ■ ko2  (1+2  (Mq  - pd  - qd))  - kx2 
C2  * 2ko2  3 


The  wave  equations  are  further  transformed  by  the  following  change 


of  variables 


5<  « - Cl1/3  z 

Ci2/3 


(■>  - + C2 V3  z 

C22/3 


■ 
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Thus  substituting  for  £<.  in  the  wave  equations,  we  have  a well 

> 

known  Stokes  equation 


(S  • •).  V 

> > 


5<d 
5 Id 

5 2d 


5ld  » $>d  * 52d 


Bl 

Ci2/3 

B2 

C22/3 


- CiV3  d 
+ C2^/3  d 


5<  * Cld  ~ Cx1/3  (z-d) 

* 52d  + C21/3  (z-d) 

and  since 

(Bi  - Cid)  * (B2  + C2  d) 
cl2/3  5ld  * C22/3  S2d 


Finally, 


®1 


510 


Cx2/3 
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1 

The  Stokes  equation  now  has  to  be  solved  with  the  help  of  the 
following  boundary  conditions: 

d Ey  dEy 

Ey  < 3 Ey  > t | 3 | Z = d 

Y dz  < dz  > 

For  perfectly  conducting  earth 


Ey  a 0 z » 0 

Ey  for  z»d  is  an  outgoing  wave.  This  condition  represents  the 
Sommerfeld  radiation  condition.  The  Stokes  equation  is  well 
known  and  it  has  solutions  in  terras  of  Airy  Integral  Functions 
Ai  (£)  and  Bi  (5)  and  they  are  given  by 

Ey<  3 a Ai  (-5<)  + b Bi  (-?<) 

Ey>  3 c (Ai  (-£>)  + j Bi  (— c> ) ) 

Here,  a,  b and  c are  the  coefficients  to  be  determined  by  applying 
the  boundary  conditions.  Note  that  up  until  now,  we  have  not 
introduced  any  modal  designate. 


Application  of  the  boundary  conditions  lead  to  the 
following  transcendental  equation 
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! 


^Ai  (- 


Ai(-Cl0) 

Cld) * Bi  (-5 Id) 


Bi(-Cio) 


5 2d)  + j Bi'  (-?2d) 


■} 


P 

- V 


1/3 


^Ai'  (- 


Ai  (-£2d)  + j Bi  (— C2d> 
Ai(-ClO) 


5ld>  - 


Bi  (-ClO > 


Bi*  (-Ud> 


) 


Here  the  prime  refers  to  differentiation  with  respect  to  S,  and 


d d d£  d 

a;  !<  " Ti  ■ 37  !<  * ' Cll/3  T?  '< 


d 1 d d?  1 , ,,  d - 

37  '»  * ai  ' 37  '>  ' + C2 v l> 


<3?  dz 

The  solutions  of  the  transcendental  equations  are  the  eigenvalues 
and  for  the  fixed  frequency,  the  wave  number  kxi  Is  related  to 
the  ith  eigenvalue  5idi  by 


kxi^  * ko^  (1+2  Mmin)  - Sidi  . Ci^/^ 


, P - /0  5ldi 

x i » (ko  (1  + Mmin)  “ ^ ^ ' — “ 


- ) 


and 


2 kQ  1/3 
( ) 


] 


I 

! 


I 

| 

1 

' 


Thus,  the  solutions  of  the  transcendental  equation  will  enable  one 
to  determine  the  height  gain  functions  and  the  eigenwave  numbers 
along  the  x axis.  Thus,  the  field  corresponding  to  the  wave  numbers 

■i 

are  completely  determined  with  one  exception,  namely,  the  amplitude 
coefficients  which  are  determined  from  the  specification  of  exci- 
tation or  antenna  characteristics. 


4.1.2  Trilinear  Refractive  Index  Profiles 

The  selective  meteorological  data  from  Canton  Island 
(Tables  2-1  through  2-12)  has  revealed  that  although,  because  of  the 
nature  of  the  measurements,  majority  of  the  cases  were  representable 
as  having  bilinear  refractive  index  profiles;  there  were  occasions 
when  the  profiles  were  more  representable  as  trilinear  and  this  in 
particular  when  the  ducts  were  elevated.  This  necessitated  devel- 
oping an  analytical  propagation  model  for  the  trilinear  refractive 
index  profile.  In  principal,  the  approach  is  identical  to  the  one 
developed  for  bilinear  case  and  it  entails  additional  boundary  con- 
ditions representing  matching  of  the  transverse  impedances  at 
the  intermediate  regions  of  M-profile. 

In  order  to  display  the  basic  features  of  the  trilinear 
refractive  index  profile,  we  have  reproduced  in  Figs.  (4-2,  3,  4) 
some  well  known  refractive  index  profiles.  The  Fig.  (4-2)  is 
obviously  a surface  duct  with  the  bilinear  M-profile,  Fig.  (4-3) 
and  (4-4)  are  in  general  categorized  as  elevated  ducts  with 
trilinear  M-profiles. 


M-Unita 


Si 


Fig.  (4-2)  - Surface  Duct  - Bilinear  Profile 


Figure  (4-3)  - Elevated  Surface  Duct  - Trilinear  Profile 
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In  all  the  three  cases,  the  duct  width  is  defined  by  d.  In  Figs. 
(4-3)  and  (4-4)  note  the  difference  in  the  elevated  surface  duct 
and  the  elevated  duct.  It  is  that  at  the  duct  apex  the  Mn,in 
is  less  than  the  M at  the  surface  for  the  elevated  surface  duct 

case,  and  Mmjn  i3  greater  than  M at  the  surface  for  the  elevated 
duct.  It  is  also  to  be  recognized  that  in  practice,  the  slope 
of  the  M-profile  will  not  be  that  sharp,  however,  the  measurement 
techniques  perhaps  would  not  provide  closely  spaced  points  to 
round  off  the  M curves. 
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4.1.3  Propagation  through  Ducts  with  Trilinear  Profiles 
The  scalar  wave  equation  for  Ey  for  the  trilinear 
refractive  index  profiles  is  given  by: 

d2 

J 2 + k02  [l+2M(zQ  - kx2  J Ey  * 0 


M(z)  Mo+qj^z 

z<di 

-*■  Mo+qidi-p  (d2~di) +p  (d2~z) 

di<z<d2 

-►  M0+q1di-p(d2-di) +q2(z-d2) 

z>d2 

or 

1 

a2  > 

_2  + b1+c12)  «*-« 

z<di 

1 

4 

["h + 82  - C2Z) Ey  ■ 0 

di<lz<d2 

and  | 

d2  ) 

+ B3  + C3Z  S Ey  * 0 

c dz2  J 

z>d2 

Here  C\j  B2, 

C2?  and  B3,  C3  are  given  by 

B1  - k02(l+2M0)-kx2 

Cl  - 2k02qi 

b2  * k02( 1+2 (Mo+qidi+pdi) )-kx2 

C2  * 2k02p 

b3  - k02(l+2(Mo+qidi-p(d2-di)-q2d2) )-kx2 

C3  * 2k02q2 

Let 

5<1  ' £3>j  +Cll/3j 

d5 

< * +C].l/3 
dz 

5<12>  - - C2l/3z 

C22/3 

d£ 

— | <12>  3 -C21/3 
dz 

b3 

£>2  ■ + C3^/3z 

C32/3 

« W3 

— |>  - C31/3 
dz 

These  substitutions  in  the  wave  equations  lead  to 


d_2 

d^2i 

2 

3 


- 0 


This  again  is  the  well  known  Spokes  equation  and 
has  solutions  in  terms  of  Airy  Integral  functions. 

Analogous  to  the  bilinear  case,  the  field  description 
in  the  three  different  regions  is  as  follows: 


Ai(-U0)  . 

Ey  a ai{Ai(-£<l)  - Bi(-5<l)} 

1 Bi(-ao) 

Ey2  - a2{Ai(-5<i2>)+a3Bi(-5<i2>)} 

Ey 3 * a4{ Ai (~C>2) +jBi (~C>2> } 

The  E field  has  to  satisfy  the  following  continuity  relations: 

Ey  ■ Ey 
Y1  y2 

dEy  dEy 

yl  » y2 

dz  dz 

Ey 2 " Ey3 

dEy2  ^By^ 

dz  dz 


These  boundary  conditions  enable  us  to  eliminate  the  different 
amplitude  coefficients. 


From  now  on  we  are  specializing  to  the  case  of  q 


q^  = q.  Consider  first  the  continuity  relations  for  field 


and  its  derivative  at  z = d^  which  lead  to  the  following  equations 


, Ai(-eio)  . 

al  {Ai(-Sldl) Bj[  (-Cldl)  ' 


Bi(>^10) 
a2  { Ai ( — C 2dl ) + 33Bi(-?2dl)^ 


and 


f , Ai(-Cio)  , , 

al  {Ai'(-Cldl) — • B±  (-€idl) > 


Bi(-5l0> 


-(£) 

q 


1/3 


*2  { Ai ' (~C2dl)  + a BiV^dl^ 


or 


Let  the  L.H.S. 
then 


= P/Q. 


p P, 


1/3 


Ai  (-C  2dl)  • Ai  ' (-Z  2dl^ 

Q q 


LBi(“^2dl)  + 


1/3 


Bi' (~^2dl) 


•^rr 
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Likewise,  the  continuity  relations  at  z - 62  lead  to  the  following 

a4  Ui(-C>d2)+jBi(-£>d2>  } " 

32  {Ai(-^2d2)+33Bi(~C2d2)l 

and  a4  { Ai ' (-5>d2) +jBi ' (~^>d2^  3 

p 1/3 

-a2  (“)  • {Ai' (-Ud2)+a3Bi*  (-^2d2)  ^ 

q 


Thus  the  transcendental  equation  is  given  by 

( Ai(-5>d2)+jBi(-5>d2)  [ . 

\ Ai' (-5>d2) +jBi ' (-C>d2) J 

/ Ai(-^2d2)+a3Bi(-^2d2) 

Ai'  (-52d2)+a3Bi'  (-52d2) 


where 


q ”1/3 
<?> 


a3 


\ Ai (~52dl)  + 

t 

{ 


p p I/3  ) 

-(-)  • Ai,(-C2dl)|L 

Q q 


/ 


P p 1/3 

Bi(-C2dl)+-(-)  Bi'(-52dl) 

Q q 


] 


4 .2  The  Solutions  of  the  Transcendantal  Equation  - The 
Bilinear  Refractive  Index  Profile 

The  transcendental  equation  which  is  arrived  at  by 
satisfying  the  boundary  conditions  at  z=d,  physically  represents 
the  matching  of  the  transverse  impedances  at  z*d. 

Zz<  * 8z> 

or  equivalently  $(2)  = (B2<-8Z>)  * 0 


E E 

-at*  JL 


Here 


and  Bz< , 8z>  are  the  transverse  impedances  for  z<d  and  z>d, 
respectively. 

Because  of  the  nature  of  the  transcendental  equation, 
a closed  form  solution  does  not  exist  and  an  iterative,  numerical 
procedure  to  locate  the  zeros  of  the  function  4>  ( B ) has  to  be 
adopted.  Three  different  approaches  have  been  explored,  and 


they  are: 


a.  The  Newton-Raphson  Method 

b.  The  Secant  Method 

c.  The  Zero  Crossing  Method 


4.2.1  The  Newton-Raphson  Technique 

This  technique,  essentially,  provides  a procedure 
to  obtain  an  approximate  solution  of  the  transcendental  equation 
which,  after  several  iterations,  approaches  the  exact  solution 


and  it  states 


HI  ■ 80  ~ 


♦ <«0> 


♦'(20) 

Here,  80  is  the  first  trial  solution  of  the  transcendental  equation 
$(Z)  * 0.  ♦'(00)  is  the  derivative  of  41  with  respect  to  80  and 

is  the  next  better  solution  to  be  used  for  second  iteration 
and  so  on.  Thus,  the  ith  solution  is  obtained  from  (i-l)th  and  is 


Bi  = Bi-l  - 


♦ (Bi-l) 
♦*  (Bi-l) 
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If  the  transcendental  equation  has  n solutions  then  this  procedure 
has  to  be  repeated  n times. 

The  Newton-Raphson  iterative  technique  to  locate 
the  zeros  of  the  transcendental  equation  is  effective  only  when 
the  trial  solution  is  reasonably  close  to  the  actual  solution. 

The  Newton-Raphson  procedure  thus  has  severe  limitations. 

4.2.2  The  Secant  Method 

This  iterative  approach  is  analogous  to  the  Newton- 
Raphson  approach  and  the  iterative  solutions  are  given  by 

(Bi-S5i-l)4>(Hi-l) 

*i  - Bi-1  

$(Bi)-<M&i-l) 

Note  that  in  the  limit  when  Bi  approaches  the  exact  solution 
8i_l,  then  this  expression  is  identical  to  the  Newton-Raphson 
expression.  Alternatively,  if  Bi-Bi-i<<<5  and  if  <}>  (« i-1) -<*>  {■)<<£ 
then  the  technique  is  identical  with  the  Newton-Raphson  technique. 
Here  6 and  e are  very  small  quantities. 

4.2.3  The  Zero-Crossing  Technique 

Another  technique  which  is  applicable  to  determine 
the  solutions  of  the  transcendental  equation  is  the  zero  crossing 
of  $ (B)  , when  <MB)  is  plotted  as  a function  of  3.  This  procedure 

is  extremely  effective  for  the  trapped  modes  or  when  the  Sldi 
are  pure  real,  positive  or  negative  numbers. 

In  any  case,  for  all  three  cases,  the  computer  programs 
for  the  Airy  Integral  Functions  and  the  Hankel  Functions  are 
required. 
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4 . 3 Numerical  Analysis  For  Bilinear  M-Profile 

It  was  decided  that  four  meteorological  profiles  will  be 
chosen  in  such  a fashion  that  the  two  extremes  of  the  character istic 
parameters  such  as  the  duct  height  and  duct  intensity  shall  be 
chosen. 

For  the  3rd  of  January,  the  intensity  of  the  duct  was 
close  to  the  maximum;  for  the  4th  of  December,  the  duct  intensity  was 
minimum.  For  the  4th  of  September  the  duct  height  was  close  to  the 
largest  while  January  3rd  it  was  close  to  the  smallest.  For  the  5th 
of  April  the  duct  was  moderately  intense.  For  all  four  days,  the 
profile  was  describable  as  bilinear.  Finally,  in  order  to  extract 
the  frequency  dependence,  three  wavelengths  were  chosen,  and  they 
were  10,  20  and  50  cm. 

Based  on  the  expression  for  the  maximum  wavelengths  that 
are  trapped  by  ducts  of  height  d,  the  following  observation  is  made: 

Table  (4-1)  - Trapped  Wavelengths 


Date 

d (meters) 

\ max  (cm) 

3rd  Jan. 

42" 

22.65 

5th  Apr. 

85 

65.20 

4th  Sept. 

94 

75.82 

4th  Dec. 

75 

54.04 

a.  For  10  and  20  cm  wavelengths  the  field  is  trapped  for 
all  the  four  cases. 

b.  For  50  cm  wavelengths  the  field  is  not  trapped  for 
Jan.  3rd  case,  however,  it  is  trapped  for  the 
remaining  three  days. 
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4.3.1  The  Zero-Cross inq  Technique  For  Trapped  Modes 


For  the  four  days  of  M-profile  data  and  for  three 
wavelengths  10,  20  and  50  cm,  the  transcendental  function  $ (B-)  is 
plotted  for  various  real  values  of  a varying  from  -5  to  +5.  These 
plots  are  displayed  in  Figs.  (4-5  to  4-16) . The  zero-crossing 
values  of  B,  or  equivalently  £ld»  ^°r  which  4>  ( B)  ■ 0 are  the 
solutions  of  the  transcendental  equation.  Consider  now  the  January 
3rd  profile  and  X ■ 10  cm  case  (Fig.  4-5).  The  zero-crossings  take 
place  for  the  following  values: 


5ld=-1.25,  0.68,  2.8,  5.25 

First  observe  that  the  zero-crossing  at  -1.67  is  not 
a true  crossing  as  is  evident  from  the  computations  for  finer 
spacings  of  the  £id  values.  The  ♦ (B)  will  continue  to  decrease 
approaching  - •<>  for  tid  values  in  the  vicinity  of  -1.65.  Likewise 
f°r  £ld-»-1.68,  the  ♦ (B)  will  increase  approach  +•;  hence  no 
zero-crossing  for  £id  * -1.67  is  possible. 

The  first  £id  for  which  ♦ (»)  - 0 is  -1.25  and  this 
negative  root  needs  to  be  interpreted.  Consider  the  field  Ey 
associated  with  this  root  and  observe  how  it  behaves  as  the  height 
increases  from  z«0  to  z«d  and  z>>d. 


Fig.  (*!-6)  Plots  of  "jfCEHs  a Function  of  ? 


Fig.  (4-7)  Plots  of^GD/as  a Function  of  2 
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Fig.  (4-14)  Plots  of  fO)  as  a Function  of  2 


Fig.  (4-15)  Plots  of  J(2)  as  a Function  of  2* 
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The  height  gain  function  or  equivalently  the  Ey  is 


given  by 


Ai(-5io) 

Ey<  = Ai  (— C<)  - ■■  ~ » Bi(-C<) 


z<d 


Bi(-ClO) 

- Ai(-5>)+jBi(-C>)  z>d 

The  C<,  £>,  have  been  defined  in  Section  4.1.1.  For  the 
Cia  value,  the  S<  range  in  values  from  +5«9  to  -1*25  and  5>  approaches 
zero  from  -3*57;  and  for  height  z = 77.8  m it  becomes  zero. 

For  greater  heights  it  is  an  increasing  positive  real  number. 

For  these  parameters,  the  field  intensity  Ey  will  behave  as  follows; 

a.  Ey  oscillates  up  to  z = 35  m at  which  £<  starts 
becoming  negative. 

b.  It  increases  in  amplitude  from  z * 35  m to  the 
duct  height  z = 42  m.  It  continues  to  inrease  up  to  z * 77.8  m 
at  which  5>-*>0 

c.  Beyond  z * 77.8  m the  field  behaves  as  if  it 
is  attenuating  with  height  since  S>  is  positive  real  number. 

The  intensity  of  the  field  above  the  duct  height  is  higher  than 
in  the  duct  and  may  be  construed  as  the  dominant  leaky  mode. 

However,  this  fact  alone  will  not  yield  the  information  as  to 
the  intensity  of  the  field  above  the  duct  since  the  excitation 
problem  has  to  be  worked  out  before  any  quantitative  estimates 
of  the  field  intensity  above  the  duct  can  be  made. 


4.3.2  The  Leaky  Modes 


Up  until  now  we  sought  solutions  of  the  transcendental 
equation  which  were  positive  or  negative  real  numbers.  The  deter- 
mination of  the  solutions  of  the  transcendental  equation  for  complex 
values  is  much  more  involved  and  time  consuming,  and  the  resources 
of  this  contract  do  not  permit  detailed  analysis.  Nonetheless, 
some  already  developed  approximate  expressions  are  used  to  determine 
the  complex  eigenvalues  or  solutions  characterizing  leaky  modes. 

The  expression  which  is  used  to  determine  the  leaky 
modes  is  given  by  (2) 


Am  - Ura+9+“2  Sg4  - ^ S5ra96 


(13S2-5S) 

1260 


1260 


SC2m98 


d P 

Here,  g*  — H*(kQ2q)-l/3f  s*l-s3,  s3»  — and  £m  are  the  zeros  of 
H q 


the  Hankel  function  of  second  kind  and  one  third  order  and  are 
given  by 


Ki  • -1. 169+j2.025 

2 -2.044+j3.54 

3 -2.76+j4.78 

The  expression  for  the  eigenvalue  Am^am+jBm 
for  small  values  of  g(g<l).  That  is,  the  duct  height  is  comparable 
to  the  scale  height  H and  for  these  values  of  g the  eigenvalues 
approach  the  zeroes  of  the  Hankel  function  of  the  second  kind 
and  one  third  order. 
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Recognizing  that  for  X=100  cm  the  field  will  not 
be  trapped  by  any  of  the  ducts;  X»100  cm  was  chosen  for  calculations. 
Thus  for  X=100  cm  and  for  two  days  of  duct  profile,  the  complex 
eigenvalues  are  presented  in  the  table  as  follows: 


Table  (4-2)  Eigenvalues  for  Leaky  Modes 


Jan. 

Dec. 

X»100  cm 

g * 0.545 

g * 0.6646 

Al 

-0.574+j2.02 

-0.48+j2.021 

2 

-1. 448+j3. 53 

-1. 354+j3. 53 

3 

-2.163+j4.78 

-2.07+j4.77 

I 

Use  of  the  same  equation  for  X»50  cm  leads  to  results 
which  do  not  appear  to  be  valid.  For  X=50  cm,  g = 1.378  for 
Jan.  3rd  duct  and  g=1.476  for  Dec.  4th  duct  and  thus  the  expression 
for  g > 1 does  not  seem  to  hold. 

4 . 4 The  Height  Gain  Functions  for  the  Trapped  Modes 

As  discussed  before,  the  variation  of  the  field  with 
height,  also  termed  the  height  gain  function,  will  reveal  whether 
or  not  there  exists  field  outside  the  duct.  For  this  reason, 
height  gain  functions  for  individual  modes  have  been  computed 
and  they  are  shown  in  Figs.  (4-17  through  4-53). 

The  height  gain  functions  are  the  plots  of  modal 
E field  i.e.,  they  are  the  plots  of  the  following  expressions: 

Ai(-ClO) 

Ey<  * Ai(-£<)  - . Bi(~£<)  z<d 

Bi (-S10) 


Ey>  * Ai(-e>)+jBi(-?>) 


z>d 
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Fig.  (4-23)  Height-Gain  Functions  - Janaary 
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Fig.  (4-38)  Height-Gain  Functions  - September 
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The  C<,  C>  and  £10  are  related  to  5i<j  where  5ids 
are  the  solutions  of  the  transcendental  equation.  The  expressions 

C<»  5>#  and  5^0  are  given  by 

£<  * 5ld  ” Ci1/3(z-d) 

£>  * 5zd  + C21/3 (z-d) 
and  £l0  * B1/C12/3 


! 


! 


Observe  the  peculiar  behavior  in  Fig.  4-18  for  heights 
greater  than  100  meters.  This  is  due  to  the  fact  that  the  arguments 
of  the  Airy  Functions  are  large  and  the  asymptotic  expressions  for 
the  Airy  Integral  Functions  are  to  be  used.  The  computer  programs 
have  not  been  modified  for  large  arguments.  The  peculiar  behavior 
displayed  in  Fig.  (4-28)  for  small  and  large  heights  is  also  due  to 
the  same  reason  of  large  arguments  of  the  Airy  Integral  Functions. 

The  height  gain  functions  which  are  shown  in  Figs. 

(4-17  through  4-53)  are  the  behaviors  of  the  individual  modes 
and  to  determine  the  total  field  these  modes  are  to  be  added  in 
phase  and  amplitude.  But  before  that,  the  amplitude  coefficients 
or  excitation  factors  are  to  be  determined  with  the  help  of  the 
antenna  characteristics  or  equivalently  the  total  field  in  the 
yz  plane  for  x»0.  The  excitation  problem  has  not  been  addressed 
in  this  report. 
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5.  SUMMARY  / RECOMMENDATIONS  I 

5.1  Summary  ' 

In  this  phase  of  the  program,  the  electromagnetic 
propagation  through  a medium,  which  was  characteristic  of  meteo- 
rological conditions  around  Canton  Island,  was  considered  in 
depth.  The  meteorological  conditions  around  Canton  Island  were 
represented  by  bilinear  and  trilinear  refractive  index  profiles. 

For  the  propagation  analysis,  both  the  geometrical  optics 
or  ray  trajectory  approach  and  the  physical  optics  or  the  guided 
wave  approach  were  considered.  The  relative  merits  of  these 
approaches  were  highlighted.  The  ray  trajactories  for  four  dif- 
ferent days  of  meteorological  data  were  also  generated. 

The  guided  mode  approach  was  pursued  to  determine  the 
propagation  behavior  for  both  the  bilinear  and  trilinear  refractive 
index  profiles. 

The  numerical  analysis  to  determine  the  modal  solutions 
and  the  corresponding  height-gain  functions  was  carried  out  only  for 
the  bilinear  profile.  The  modal  solutions  were  obtained  by 
iterative  numerical  analysis  of  the  transcendental  equation  for 
the  trapped  modes  only. 

For  the  leaky  modes,  a previoushy  developed  approximate 
expression  was  used.  The  validity  of  the  expression  for  duct 
heights  smaller  than  the  frequency  dependent  scale  height  was 
established. 


> 
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Although  leaky  modes  are  indicative  of  the  field  above 
the  duct  height,  for  quantative  assessment  of  the  field  intensity 
in  the  shadow  zone,  the  excitation  problem  has  to  be  solved. 

This  has  not  been  addressed  in  this  phase  of  the  program. 


i 


5.2  Recommendations 

The  usefulness  of  the  guided  wave  propagation  analysis 
hinges  on  the  ability  to  generate  radar  coverage  diagrams.  The 
radar  coverage  diagrams  in  the  presence  of  different  kinds  of 
meteorological  ducts  would  enable  one  to  evaluate  the  performance 
degradation  of  ESD  sponsored  surveillance  radar  systems. 

For  the  generation  of  radar  coverage  diagrams,  the 
following  tasks  need  to  be  completed. 

1.  An  effective  analytical  technique  to  determine 
the  complex  eigenvalues  representing  leaky  modes. 

2.  Determining  the  amplitude  coefficients  of  the 
trapped  and  leaky  modes  by  addressing  the  excitation 
problems  or  by  specifying  the  antenna  characteristics. 

3.  Applying  these  techniques  to  different  meteoro- 
logical conditions  and  different  radar  systems 
parameters. 

These  tasks  can  constitute  the  second  phase  of  the 

program. 
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APPENDIX  I 

Computer  program  to  determine  the  roots 
of  the  transcendental  equation. 
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